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Generalized universality, as recently proposed, postulates a universal non-Gaussian form of the 
probability density function (PDF) of certain global observables for a wide class of highly correlated 
systems of finite volume A''. Studying the 2D XY-model, we link its validity to renormalization group 
properties. It would be valid if there were a single dimension operator, but the actual existence 
of several such operators leads to T-dependent corrections. The PDF is the Fourier transform of 
the partition function Z(g) of an auxiliary theory which differs by a dimension perturbation with 
a very small imaginary coefficient iqlN from a theory which is asymptotically free in the infrared. 
We compute the PDF from a systematic loop expansion of InZ(g). 
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Trying to understand universality properties of crit- 
ical systems near a second order phase transition was 
extremely fruitful for the development of physics. It lead 
to the development of Wilsons renormalization group 
(RG) which found applications ranging from solid 
state physics to elementary particle theory. Universal- 
ity classes in the RG-sense depend on the dimension of 
the system and on the symmetry properties of the order 
parameter. 

Recently a generalized universality has been pro- 
posed 0, which goes far beyond the known picture. It is 
therefore of great general interest to understand the con- 
ditions for its validity. Generalized universality is sup- 
posed to hold true for systems sharing the properties of 
strong correlations, finite volume and self similarity, no 
matter whether they have the same symmetries and di- 
mensions nor whether they correspond to equilibrium or 
non-equilibrium systems. Such new universality is ex- 
pressed in a non-Gaussian universal curve these systems 
share when the PDF of some global quantity of each sys- 
tem is plotted (universal fluctuations). Studying the 2D 
XY-model as an example, we link its validity to RG- 
properties of the system, viz. the existence of dimension 
perturbations of a system which is asymptotically free 
in the infrared. It would be exactly true if there were 
one single dimension operator. But the actual exis- 
tence of several dimension operators in the model leads 
to T-dependent corrections to the supposedly universal 
curve which do not go away when the volume iV i-^ oo. 



•Electronic address: Igerhard . mack@ desy. de| 
^Electronic address: gpalma@lauca.u sach.cl| 
■f Electronic address: lvergara@lauca. usach.cl^ 



If there were no dimension operator, the PDF would 
be Gaussian. 

Let us recall some previous results. The supposed uni- 
versal curve for the PDF, which we call BHP, is non- 
Gaussian in spite of the result that a naive application 
of the central limit theorem would have given; this pos- 
sibility has been attributed to the anomalous statisti- 
cal property of finite size critical systems that, whatever 
their size is, cannot be divided into mesoscopic regions 
that are statistically independent, a necessary condition 
for the central limit theorem to apply. In Ref. Q it was 
shown that a steady state in a closed turbulent flow ex- 
periment and a finite volume magnetic model, the 2D 
XY-model in a finite square lattice, had the same uni- 
versal non-Gaussian probability distribution function for 
the fluctuations of a global quantity. These are two ex- 
amples of systems, among others (see 0, O, 0, 0), 
that exhibit such data collapse. 

In it was argued that in finite-size systems that are 
in the critical regime standard scaling is a sufficient con- 
dition to exhibit data collapse. Nevertheless, as our re- 
sults show, when one computes instead of the PDF more 
sensitive quantities (from the numerical point of view), 
such as its moments (skewness and kurtosis), tempera- 
ture dependence shows up. This occurs within the range 
of temperatures of the expected data collapse. 

A functional form for the BHP universal curve was sug- 
gested in and it was also argued that this universal 
curve should be independent of both size and tempera- 
ture of the system. This result was based on an earlier 
computation of the PDF of the magnetization of the 2D 
XY-model in the harmonic (=spin-wave) approximation 
(the 2D HXY-model) 0. In this reference the PDF of 
the magnetization was carried out by resumming its mo- 
ment series up to leading order in N (the system size). 
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assuming that multiple loop diagrams can be neglected, 
and finally by performing an inverse Fourier integral. 

Nevertheless, more recently ^3 a numerical study of 
the full 2D XY-modcl was performed, and the results 
suggested a small but systematic dependence of the PDF 
on the system temperature. This numerical analysis was 
performed by using high precision Monte Carlo simula- 
tions, which amounts to control the statistical indepen- 
dence of the configurations used to compute thermody- 
namical averages, as this system is affected by the well 
known critical slowing down effect. 

In this work we explain how to compute systematically 
the PDF for the 2D XY-model within the loop expan- 
sion of an auxiliary theory with an imaginary coupling. 
It turns out that the 1-loop approximation agrees with 
the curve BHP computed in Ref. 0| . But higher loop 
corrections have a different T-dependence. They are not 
suppressed by powers of and the T-dependence of 
their contribution to supposedly universal quantities such 
as skewness and kurtosis does not disappear for large N. 
This unexpected breakdown of strict BHP universality is 
not due to the corrections to the harmonic approxima- 
tion to the Hamiltonian, but comes from the fact that 
the nonlinear order parameter (the magnetization) is a 
sum of several dimension operators which contribute at 
different order of the loop expansion. 

We will later confirm the known fact that the cor- 
rections to the spin wave approximation will not affect 
the PDF, for temperatures well below the Berezinskii- 
Kosterlitz-Thouless critical temperature Tbkt and for 
large iV, apart from a renormalization of temperature. 
Therefore the PDF for the XY and the HXY-model 
have the same large N limit except for a renomalization 
of temperature, but neither of them obeys strict BHP- 
universality. 

Because vortices are neglected in perturbation theory, 
results are only valid well below Tc 

Let us consider a classical statistical mechanical sys- 
tem which lives on a lattice with N sites x, to which 
variables (j)x are attached, and with a Boltzmann weight 
g-/3W(0)^^^ where Dcf) = Ilx d4>x and dcj)^, IS a measure 
on the space in which (ji^ takes its values. Let be a 
global quantity of the form 

X 

where '^{4>x) is some function of ipx- The magnetization 
per site in a ferromagnet is an example of such a global 
variable. 

The PDF is defined as 

P{M) = - / D4>e-P'^'^'>'H{A{(I>) - M) (1) 
Za J 

with partition function Zq — J D4>e^^^'^'^'^\ The mean 
(M) and mean square fluctuation a'^ = {{M-{M)Y) of A 
are computable from the PDF as (M) = / M P{M)dM 
and cr^ = / {M - {M)fP{M)dM, respectively. 



Inserting the Fourier representation of the 5-function, 
we have 

P{M)^ r ^e-^«(*^-<*^»Z(g), (2) 
J-oo 27r 

where 

Z{q) - j D^exp | -/37^(0) + ^ ^ $(0.) | . 

(3) 

The hypothesis of BHP universality says that the PDF, 
considered as a function 

P(M) ^ F [''-J''^) 

of {M — {M))/a, is a universal function F, the same for 
many strongly correlated critical systems, independent of 
temperature T = and of system size N, provided N 
is large enough. 

We see that the PDF is the Fourier transform of the 
"partition function" Z{q), normalized to Z{0) — 1, of an 
auxiliary theory whose Hamiltonian differs from that of 
the original theory by a perturbation with an imaginary 
coupling iq/N which is small of order 1 /volume = N~^. 
On finite systems, Z{q) will typically be a holomorphic 
function of g in a small neighborhood of zero, and Z{q) 
is a true partition function for positive imaginary q. We 
choose to extract the factor exp{iq{M)) from the parti- 
tion function so that its derivative Z'{0) = 0. 

The hypothesis of BHP-universality asserts that for the 
above mentioned systems, the "partition function" Z(q) 
is a universal non-Gaussian function of x = —aq. We 
consider the magnetization per site of the 2-dimensional 
XY-model. (j)x can be considered as real variables, and 
the Hamiltonian is 

= -j - •^J')- (4) 

<x,y> 

Without loss of generality, the coupling constant J > 
can be set to 1. The sum runs over pairs of nearest 
neighbor sites x,y on a square lattice of side length 
L = \/N . The Hamiltonian is invariant under a global 
rotation of the spins (cos sin 0x) and under transla- 
tions (j)x ^ ipx + 2Trnx, Tlx integer. By a global shift of 
the variables <f>x one can achieve X^a; '/'a; ~ ^- The mag- 
netization per site is then given by 

$(0^) = COS02;. (5) 

The density cos <j)x is a sum of dimension perturbations 
which are eigenmodes of the linearized renormalization 
group. The existence of a dimension perturbation ex- 
plains why the (/-dependent perturbation can make a con- 
tribution of order to In Z{q), in each order g", despite 
the small coefficient cx . The existence of several (in 
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fact infinitely many) dimension perturbations explains 
why there is no strict BHP-univcrsality. If the pertur- 
bation had dimension > 0, successive orders in q would 
be suppressed by negative powers of N, so that the term 
cx would dominate after the term cx g is subtracted 
from In Z{q), and the PDF would be Gaussian. 

The existence of several dimension operators is un- 
usual even in two space dimensions. But for the 2D 
XY-model it follows from features of free field theory 
because the (known) accuracy of the harmonic approxi- 
mation to the Hamiltonian implies that we consider RG- 
transformations in the vicinity of a Gaussian fix point. 

In the spin wave approximation, the Hamiltonian is 
approximated by 



X fl — 1 



(6) 



where is the finite difference derivative. This is the 
Hamiltonian of the HXY-model. 

Since /3 = 1/T appears as a factor in front of H, the 
logarithm of Z{q) can be computed as a power series in 
T, for fixed qT/N . This is called the loop expansion since 
the parameter T plays the role of ft in field theory [T^ . 
The expansion up to n loops includes terms up to order 

rpn — l 

On finite lattices, InZ(g) may be expanded in a power 
series in q. 



n>l 



The coefficients a„ also depend on system size N . By 
definition, Z{q) is normalized to 1 and obeys Z'{Q) = 0, 
hence 



1 - (M) = Tai{T) 



(8) 



and the mean square fluctuation cr^ = —-;pr^ \nZ{q = 0) 
equals 

cr^ = 2T^a2{T). (9) 

To obtain the PDF with normalized first and second mo- 
ments, one considers InZ(q) as a function of x = —qa. 
It has the form 

-iq{M - (Af)) + \nZ{q) = ix^^^—^ - ^ + ^W, 

with ^'(a;) — 0{x^), and the higher moments are ob- 
tained from it as 



M - (M) \ " 



1 2 

2 



(10) 

for n > 3. Inserting the power series expansion ((T)) of 
In Z into eq. (0) results in 



dx 

P{M) = / - — exp \ ix 



{M - (M)) 1 



n>3 



(11) 



The skewness s = {{M — {M))^)a ^ and kurtosis c 
{{M - {M))'^)a^'^ are obtained from eq.lO as 
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T 



as 



T 



"4 



(12) 
(13) 



BHP-universality would require that the coefficients 
an{T){T / aY^ of {ix)" have universal values, given by the 
BHP-formula g„(ig2)""/V2n which we recover below as 
the 1-loop approximation. In particular, these coeffi- 
cients should be T- independent, and also the skewness 
and kurtosis should be T- independent. We will compute 
the coefficients a„ to 2-loop order. 

We extract the quadratic terms 



/3H(0) = msw{4>) + Vi{4>), 



(14) 
(15) 



COS( 



= f E(l-^'^')-^2(^.) (16) 



N 2 

X X 

and define the q-dcpendent propagator 
F 



(/ + ^G)-iG, 



iq{M)+\nZ{q) = iq+^{-iqTYan{T) (7) where 



G{x) 



G(k) = ^ 



j^-l^^-ikx 
1 

K 

E sin'(fcM/2) 



G(k), 



(17) 

(18) 

(19) 
(20) 



is the propagator of a massless scalar free field on a finite 
2 dimensional lattice with the zero mode removed. 

Using the familiar device of introducing a source term 
(j, (/)), the integral defining Z{q) can be formally per- 
formed with the result 



Z{q) = - e-'"^^^'^-^'' exp (iTrln[TF] 
Zo I 2 



• exp < —V 



T 



(21) 



3=0 



with V = Vi+ V2. 

In 1-loop approximation, we obtain from Eq. (|21|l 



-iqiM ~ (M))+lnZ{q) 



-iq{M -1) + ^Tr ln[F] 



plus g-independent terms. Upon inserting the definition 
(113 of F, 

TrlnF = ^ |lnG(k) - In (l + ^G(k)) | (22) 



4 



the PDF in 1-loop approximation can be computed by 
numerical Fourier transformation according to eq. Q • A 
similar computation is also possible when the 2-loop cor- 
rection is added; in harmonic approximation it requires 



from the (j)'^ contribution to ^ cos < 

-iq{M- {M))+ In Z{q) is 



Its contribution to 



(27) 



r(0) = >F45:(i + !fG(k 



G(k) 



according to ea. l|77| below. 

This result of the PDF in 1-loop approximation repro- 
duces exactly the previous result 0], 0| for the BHP 
PDF. In section II A of it is stated It appears that 
the values of multiple loop graphs, like the first one in 
Fig. Ic), are zero in the thermodynamic limit". But in 
fact this is not true for all graphs, and as a result the 
BHP PDF is merely a lowest order approximation which 
becomes exact in the limit T i— )■ only. 

We will compute the 2-loop contribution below and 
comment on the precise relation between the diagram- 
matic expansion here and in |T]| later on. 

Following Ref. we define the quantities gn via pow- 
ers of G as 



G"(0) = N-'^TrG" = N"-'^gr, 



(23) 



Expanding F in powers of q gives 

k-l 



a 



{2harm) 



T 



-y 

8 ^ 



gm+igk- 



(28) 



This term is also present in the HXY-model. 

The second correction comes from the anharmonic 
term —(3{\7(j)x)'^/4:l in the Hamiltonian. It is absent in 
the HXY-model. It evaluates to 



|iV ([Ar(0)]2 - [AG(0)]^ 



(29) 



The g-independent subtraction cx [AG(O)]'^ restores the 
normalization In Z{0) = 0. Expanding in powers of q and 
including terms of order 1 /N results in 



(2anh) 



-y 



9k-igi 



(30) 



1=1 



gi is logarithmically divergent in the thermodynamic 
limit, gi - i^lnGiV, G = 1.8456 [ll| while the re- 
maining ones become A'^-independent for large enough N 
and in the thermodynamic limit they read 



52 = 0.38667 lo- 



gs = 0.7572 10" 



34 = 1-763 10' 



Expanding the q-dependent propagator in powers of q 



we obtain the 1-loop approximation a 
ficients a„(T), viz. 



4^nT) = |^ 



(1) 



2n 



(T) of the cocf- 



(24) 



together with the results for the expectation value (M) 
of the magnetization and its mean square fluctuation 
in 1-loop approximation 



(M)(i) = l-Tgi/2 
(a(i))2 = T'g2/2. 



(25) 
(26) 



Note that a is proportional to T, so that the factors 
(T/a)'' are T-independent in 1-loop approximation. 

Inserting the result 124(1 into ea. Hll|l we recover the 
known formula of BHP for the universal PDF. 

Numerically, the skewness s — {{M — {M))^)a~^ and 
kurtosis c = {{M — {M)Y)a^'^ in 1-loop approximation 
are obtained from Eq. as s = 0.8907 and c = 4.4150 
for I— > oo. 

In 2-loop approximation there are two corrections, 
Un = ^-a"n^°'™^ -'rOL"n°'"'^\ The first Correction comes 



We shall argue below that to leading order in 1/N this 
anharmonic correction only amounts to a renormalization 
of temperature which cannot lead to a violation of BHP- 
universality. 

The presence of the 1 /A^-corrections may explain why, 
for finite lattice sizes, the shape of the PDF of the full 2D 
XY-model displays a stronger temperature dependence 
than the spin- wave approximation, as shown in Ref. |lOj . 

In 2-loop approximation, to next to leading order in 
T, the average magnetization, mean square fluctuation, 
skewness s and kurtosis c including anharmonic correc- 
tions come out as 



(M) 



T 



:9i 



-51 



9l 



(31) 



,2=T^||l_Tgi+r (l + A.|)^, (32) 



-53 



A. 

252 



= 3 1 



52 



3/2 



l + T 



4 53 



3 

27V 



5i52 

53 



(33) 



454 

9 

52 



4T 

52 



- 252.93 + 



25153 



2 5?.94 

"-52 

52 



(34) 
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Using the values of f/„ for the appropriate volumes N, 
skewness and kurtosis for the XY-mode\ evaluate to 

s = -0.891 + 0.1319 T (N^oo) (35) 

= -0.854 + 0.1555 T (iV = 16^) (36) 

c = 4.415- 0.470 T (iV = oo) (37) 

= 4.328 -0.504 T (TV = 16^). (38) 

As stated above, terms suppressed by powers of N^^ 
come from the anharmonic corrections. 

We see that both s and c have a T-dependence which 
docs not disappear as N oo. This disproves strict 
BHP-universality. 

The values for N — oo and N — 16^ are close to each 
other, yet different enough to suggest that comparison 
with Monte Carlo data at = 16^ will yield a sensitive 
test of the accuracy of the calculation. 

In Fig. 1 we show Monte Carlo results, similar to those 
of Ref. ^3 1 for the skewness and kurtosis in the tem- 
perature range T — 0.2, ...,1.8 together with the T = 
result given by the 1-loop approximation, and the 2-loop 
results (jSnj, up to order T. The agreement is quite 
satisfactory. 




FIG. 1: Skewness and kurtosis. Full triangles and circles are 
Monte Carlo data for a lattice with A'^ = 16^ points. Full lines 
correspond to the numerical values given in Eqs. 11361381 . The 
exact 1-loop results at T = are shown as full square dots. 
Dotted and dash-dotted lines correspond to the thermody- 
namic limit {N — oo). 

In 1-loop approximation, the (iq/N) cos (j)x- 
perturbation is approximated by [iq/N)[l — 0^/2). 
The 2-loop correction takes the 0!|-term in the Taylor 
expansion of cos 0a; into account. The different T and 
cc-dependence of the 2-loop correction to In Z{q = —x/a) 
can therefore be attributed to the fact that it comes 
from a different relevant operator of dimension 0. 

We will now argue that the effect of the corrections 
to the spin wave approximation on the PDF is merely 
a renormalization of temperature, up to corrections of 
order N~^, viz. 

Th^r(i + r/2 + ...) 



This comes from the fact that the theory with Hamilto- 
nian Ti is asymptotically free in the infrared. 

The 1-loop approximation is not affected by the renor- 
malization of temperature, because it is T- independent, 
and the 2-loop correction is also unaffected to order T, 
because the correction would be of order . Thus, to 
order T, the skewness and kurtosis in the XY and the 
HXY models differ only by corrections proportional to 
1/A^. 

Because the coefficient of the g-dependent perturba- 
tion in Z[q) is of order iV~^, the g-dependent part of 
In Z{q) is entirely determined by the infrared behavior 
of the theory, up to corrections of order N~^. There- 
fore it could be computed from an effective low energy 
theory. In principle its Hamiltonian would have to be 
computed by a real space renormalization group (RG) 
calculation [l| . Since Vi is a sum of irrelevant perturba- 
tions of massless free field theory, their coefficients get 
suppressed when the length scale gets enlarged, but they 
can give rise to a renormalization of the coefficient 
of the marginal |V(^a;p term. This is the aforementioned 
renormalization of temperature. In principle there could 
also be a renormalization of the coefficients of the rele- 
vant perturbation, i.e. of the coefficient q of (j)^ in 1-loop 
approximation. But a multiplicative renormalization of 
q drops out when the second moment of the PDF is nor- 
malized. New terms in the perturbed effective action like 

cos 20a; can appear in higher orders, but they are too 
small to be detectable. 

Note, however, that the perturbative treatment of the 
corrections to the spin wave approximation does not take 
vortices into account. Its results are therefore only valid 
well below Tc. 

To verify the RG-argument, we inspect the anharmonic 
corrections (150]) . Addition of Tgk/4: to 
corrects the 1-loop result by the substitution 

gk^ 9kil + kT/2) (39) 

to first order in T. Clearly, this has no effect on gk<J^^ 
to order T, and therefore the normalized PDF remains 
unchanged. 

Next we discuss the precise relation between our loop 
expansion and the diagrammatic expansion of Ref. |ll| , 
and hyperscaling. 

One may envisage generalization of the kinetic term 

f EE I^M'/'.P = -f ('^, A0) = ^(0,G-V) 

by admitting arbitrary positive G, i.e. arbitrary prop- 
agator G{x,y). [In the translation invariant case, the 
propagator is considered as a function G{x — ?/) of a sin- 
gle variable]. Diagrammatic expressions exhibit results 
in their dependence on G, as sums of products of propa- 
gators G{x, y). Therefore they have a universal meaning, 
and the diagrammatic expansion for the same quantity 
must agree, no matter how the computation is organized. 
Our way of organizing the expansion is more economical 
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in that it avoids having to go through the moments; they 
are here summed from the start and Z(q) is their gener- 
ating function. We obtain diagrammatic expansions for 
(connected parts {M^Y of) {KP) by expanding our q- 
dependent propagator T and Tr hi F in powers of g, cf. 
eas. H17(l . H22(l and 123|l . To recover the diagrammatic ex- 
pansion of Ref. 11] for {Mp)/{M)p one must divide by 
{M)P. As shown in the Appendix, the loop expansion can 
be reorganized so that (M) appears as a factor multiply- 
ing q. This implies that the diagrams for (MP) / {M)p are 
precisely the diagrams for (MP) which contain no tad- 
pole. A tadpole is a line in the graph whose source and 
target agree; tadpoles produce factors TG{0)/2 = Tgi/2. 
Figure 2 shows the connected 2-loop graphs without tad- 
poles which contribute for p = 3, 4, in spin wave approx- 
imation. 

b) 

FIG. 2: Connected 2-loop graphs in spin wave approximation, 
a) contribution to i^{M^Y/'i\ = (— j(jr)^Q3, b) contribution 

to i^{M*Y/^\ = (-igr)V. 

For the sake of comparison let us evaluate these dia- 
grams. 

We remember that in a scalar field theory with inter- 
action V = Ad)^. vertices have n incident lines and con- 
tribute a factor —A. The —sign comes because the Boltz- 
mannian involves . In the expansion of -^cos^^:! 
even powers 0" appear with factors 1/n! and alternat- 
ing sign. The Feynman rules for our theory in harmonic 
approximation say that iq{M) + In Z{q) is the sum of all 
connected vacuum diagrams, and in these diagrams there 
is 

• a factor ]^(— )"^^ for every vertex with an even 
number n of incident lines 

• a factor TG{x — y) for a line from a vertex at x to 
a vertex at y 

• a sum over x for every vertex 

• a symmetry factor where g is the number of 
elements of the symmetry group of the graph. 

Since every vertex contributes a factor q, the vac- 
uum diagrams with k vertices sum to the summand 
{—iqTYak{T) in the expansion Q. Contributions with 
tadpoles can be eliminated by setting 171 = 0. The sym- 
metry factor for the graph (a) is | , and the symmetry fac- 
tor for the graph (b) is i. Using that J^y ~ y)G{y — 



x) = 0^(0) = Ng2 etc, the diagrams evaluate to results 
in agreement with ea. H28|) . viz. 



3 ^ 2harrn 



{-iqTYai 



-iqTfaf""'' 



i^qfY^l + o{gY (40) 
~i^qf^g293 + o{gY ■ (41) 



We see that the contribution of the diagrams shown in 
Fig. 13 does not vanish in the thermodynamic limit. 

The result that (Af ) can be factored out confirms that 
hyperscaling remains true to all orders in the loop expan- 
sion. Hyperscaling asserts that (M)^/(t^ has a thermo- 
dynamic limit. Within the frame work of a loop expan- 
sion, this must be true order by order in T. This is not 
trivially true because gi has no thermodynamical limit. 
Hyperscaling also remains true in the presence of anhar- 
monic corrections because these merely renormalize the 
temperature, as we saw. Our 2-loop results for a and 
(M) which include anharmonic corrections give 



T{M) (1 + -T + 0{T^) 



(42) 



for large N, which should be compared with eq. (19) of 
Ref. O. 
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Appendix: Normal ordered procedure 

We introduce the normalized Gausssian measure with 
covariance TG 

rfMTG(</') = ^^exp|-^^|V^^,p|Dc/, (43) 

It has the characteristic function 

d/iTG(0)e'<^'''^> ^e-^^^-^^K (44) 

If A{(f>) is a function of which is integrable with respect 
to the Gaussian measure, it can be written in normal 
ordered form 14] . The normal ordering operation is typ- 
ically indicated by ::, but it depends on the covariance of 
the Gaussian measure. 



Ai^) = : 5(0) :, 



In particular 



-ia"TG(0) . ^tai 



(45) 
(46) 

(47) 
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It follows that 



Therefore 



COS( 



= e-^™(") :cos, 



diJ.Tci'f) : cos^j. := 1 . 



(48) 
(49) 



Under general conditions, i? is a holomorphic function 
of (j) which can be expanded into a power series. In par- 
ticular 



COS( 



1 



1 

'4! 



(50) 



It is well known that the normal products : 0^ : 
are always eigenmodes of the linearized renormalization 
group at the Gaussian fix point determined by TG. In 
two dimensions, with fix point Hamiltonian pjisw — 
i/?^ |V^(/)p they all have dimension 0. 

Consider now the partition function 



Ziq) = Z((M)q) 



(55) 



with a new partition function 



Z{q) = J c?MTG(0)exp 1^ ^ : cos :| , 



(56) 



and 



(57) 



Z(q) = e^^<'''^Z(q). (51) 
It follows from cq.(|21l that the moments 

{MP) = (^-^^J Ziq)U=o. (52) 

We restrict our attention to the spin wave approxima- 
tion, replacing H by Hsw- Inserting ea. H48|l into the 
definition Q of Z{q), we obtain 

Ziq)^ J dMTG(0)exp||e-?«W^:cos0, :| 

(53) 

From this and eq. 149|) we recover the known result for 
the magnetization in spin wave approximation, 

(M) = e-i'^^^^ (54) 



The partition function Z{q) is exactly like Z{q) except 
for the normal ordering of cos ■ In a diagrammatic 
expansion, normal ordering eliminates tadpoles. The di- 
agrams which contribute to (M^) / {M)p are therefore ex- 
actly the diagrams without tadpoles which contribute to 
(MP). 

It is more convenient to consider truncated expectation 
values 



(MP)7(M>P= lnZ(g)|,=o. (58) 



They have expansions in connected diagrams. The n- 
th order in the loop expansion for \nZ sums connected 
diagrams with n loops. In 1-loop approximation only 
the term 1 — i : 0^ : in the expansion (|50|l of : cos 4>x ■ 
contributes, but the jf : <j)'^ :-term begins to contribute 
in 2-loop order, and so on. 
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